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Using atomistic simulations we determine the roughness and the thermal properties of a suspended 
graphane sheet. As compared to graphene we found that hydrogenated graphene has: 1) a larger 
thermal contraction, 2) the roughness exponent at room temperature is smaller, i.e. ~ 1.0 versus 
~ 1.2 for graphene, 3) the wave lengths of the induced ripples in graphane cover a wide range 
corresponding to length scales in the range (30-125) A at room temperature, and 4) the heat capacity 
of graphane is estimated to be 29.32±0.23 J/molK which is 14.8% larger than the one for graphene, 
i.e. 24.98±0.14 J/molK. Above 1500 K we found that graphane buckles when its edges are supported 
in the x — y plane. 



I. INTRODUCTION 

Graphane (GA), a two dimensional covalently bonded 
hydrocarbon, was first predicted from ab-initio calcula- 
tion by Sluiter et al [l[ and recently rediscovered 2]. In 
a recent experiment Ellas et al Q demonstrated the fab- 
rication of GA from a graphene (GE) membrane through 
hydrogenation which was found to be reversible. Den- 
sity functional theory and molecular dynamics simula- 
tions employing different force fields were carried out to 
study the structural and electronic properties of both GA 
and GE 

Hydrogen (H) atoms are chemically bound to the car- 
bon (C) atoms on alternating sides of the membrane 
(chair like conformer) which causes a local buckling of the 
membrane. Such deformations for small membrane sizes 
(< Inm) have been recently reported @. In the early 
stages of the hydrogenation process, membrane shrink- 
ing and extensive membrane corrugations occur due to 
the formation of a significant percentage of uncorrelated 
H frustrated domains Q. 

The morphology of perfect GA, and its comparison 
with a perfect GE membrane, has not yet been investi- 
gated for large samples. Ab-initio calculations are re- 
stricted to small unit cells and therefore we will use 
atomistic simulations to show the main differences be- 
tween the morphology of a large sample of GA and GE. 
Our atomistic simulations are based on the second gen- 
eration of Brenner potential Q (REBO) which includes 
the interaction with third nearest neighbors via a bond- 
order term that is associated with the dihedral angles. 
Therefore, such a potential is suitable for various atom- 
istic simulation purposes including the calculations the 
lattice thermal properties (phonon dispersion Q, elas- 
tic moduli thermal conductivity [l3|, etc) of carbon 
nanotubes, graphene and hydrocarbons. Nevertheless, 
there are a few shortcomings in REBO potential which 
are important when modeling processes involving ener- 
getic atomic collisions (because both Morse-type terms 
approach finite values at very small atomic distance), and 
the neglect of a separate tt bond contribution which leads 



to the overbinding of radicals our study is insensitive for 
these shortcomings). Of course, the electronic proper- 
ties of GE and GA and the thermodynamical properties 
at low temperatures (because of the e.g. quantum zero 
energy) are beyond REBO's ability. In those situations 
ab-initio calculations, particularly those based on den- 
sity functional theory (DFT), are extremely useful. Ab- 
initio molecular dynamics simulation (e.g. Car-Parrinello 
molecular dynamics [Tlj ) is a highly appreciated theory 
which eliminates the force-fields based restrictions, but 
its main disadvantage is that it is only applicable for 
small size systems (typically N < 500 atoms). 

Here we study relatively large systems and con- 
sider thermal and structural properties of GE and GA 
above 50 K and show that for a GA sample with size 
183 Ax 185 A containing 22400 atoms, the roughness, the 
induced ripples structure and the total energies are very 
different from a GE membrane. It has been shown that 
the ripples in GE strongly affect its thermo-electronic 
properties [l^l- Temperature effects are studied and we 
found that the thermal contraction of a suspended GA 
is larger than the one for GE. Moreover the calculated 
roughness exponent indicates that GA is more rough than 
GE even at room temperature. From our simulations, we 
predict that the ripples in GA affect its thermo-electronic 
properties much more than for GE. 

This paper is organized as follows. In Sec. II we in- 
troduce the atomistic model and the simulation method. 
Section HI contains our main results for both graphene 
and graphane. Results for thermal contraction, rough- 
ness and heat capacity for different temperatures are pre- 
sented and compared to available experimental results. 
In Sec. IV we conclude the paper. 



II. METHOD AND MODEL 

Classical atomistic molecular dynamics simulation 
(MD) is employed to simulate large flakes of GE and GA 
at various temperatures. The second generation of Bren- 
ner's bond-order potential is employed which is able to 
describe covalent sp^ bond breaking and the formation of 



2 




FIG. 1: (Color online) Two snap-shots of suspended graphane (GA) and graphene (GA) at 700 K. The inset is a zoomed region 
which shows a regular alignment of C-H bonds. 



associated changes in atomic hybridization within a clas- 
sical potential 7]. The Brenner potential (REBO) terms 
were taken as Ep = J^i J^jyii^^itj)^ BijV^{rij], where 
Ep is the average binding energy, and and are 
the repulsive and attractive term, respectively, where r.^ 
is the distance between atoms i and j. Bij is called the 
bond order factor which includes all many-body effects. 
Bij depends on the local environment of the bond i.e. 
the bond and torsional angles, the bond lengths and the 
atomic coordination in the vicinity of the bond. This 
feature allows the Brenner potential to predict correctly 
the configurations and energies for many different hydro- 
carbon structures and the correct hybridization. 

Initially the coordinates of all carbon atoms in GA 
(GE) are put in a flat surface of a honey comb lattice 
with nearest neighbor distance equal to Oq = 0.153 nm 
(0.142 nm). The hydrogen atoms are alternatively put 
on top and bottom of the carbons, i.e. (CiJ)„ with 
n — 11200 (chair like model for GA 0, see the inset 
in the left panel of Fig. [1|). To make sure that the sec- 
ond generation of the Brenner potential gives the correct 
force field for our thermo-mechanical studies of the GA 
sheet, we performed extra simulations where we started 
with a non optimum bond length of gq = 0.142 nm. Al- 
ready after only a few time steps the correct optimum GA 
lattice spacing, i.e. oq ~0.153 nm, was found confirming 
the ability of the used potential for GA simulations. 

The considered systems are square sheets with dimen- 
sions 183 Ax 185 A for GA and 170 Ax 170 A for GE, in 
X and y- directions where we considered both armchair 
and zigzag edges. The number of atoms in GA (GE) is 
n — 22400 {n/2). We simulated the system at non-zero 
temperatures (above 50 K) by employing a Nos'e-Hoover 
thermostat. In order to mimic the experimental set-up 
we prevent motion along the z-direction at the ends of 
the system along x. This is realized in practice by fixing a 
row of atoms at both longitudinal sides. This supported 
boundary condition prevents the drift (of the ends of the 



sheet) in the z-direction. Motion in the x — y plane is 
allowed which allows the system to relax and to exhibit 
lattice contraction/expansion. 

Before starting the sampling, we let the system to find 
its true size and temperature during the first 5 ps. During 
this process GA and GE change their longitudinal length 
so that we always observe lattice contraction implying 
that GA and GE shrink through surface corrugation (i.e. 
creating ripples) [l^ . Figure [1] shows two snap-shots of 
GA and GE relaxed at 700 K. The C-H bonds (except 
at the boundaries) are almost parallel which is a conse- 
quence of lateral H-H repulsion. 



III. RESULTS AND DISCUSSION 

A. Equilibrium length and thermal contraction 

After equilibrating, we start to record the longitudi- 
nal size (in x — y plane along x-axis) of the system for 
2 ps. Averaging over the lengths during the sampling 
process gives the average length of the system, i.e. (L), 
at the given temperature. Figure [H^a) shows the ther- 
mal contraction coefficient, i.e. = x^^j versus tem- 
perature, where Lq ~ 183 A (io ~ 170 A) is the ini- 
tial length of GA (GE) (for a ffat sheet at T^O K) and 
dL = (L) — Lq. As we see 7 increases with temperature 
and is always negative in the studied temperature range. 
Surprisingly GA has a larger contraction with tempera- 
ture than GE. These curves are fitted to a -I- (3/T where 
a=22±5xlO-6K-i, 4±lxlO-^K-i and /3=-877.10±29 
and -298. 8±6 for GA and GE, respectively (solid curves 
in Fig. 2(a)). 

At present a few studies exist on the negative thermal 
expansion coefficient (TEC) of GE, but none are available 
for GA. The actual values of the TECs for GE differ from 
one study to the other depending on e.g. the size and 
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the boundary conditions [T^l- Mounet and Marzari [T^ 
presented a study of the thermodynamical and struc- 
tural properties of carbon based structures (graphite, di- 
amond and GE) using GGA-DFT calculations. In or- 
der to calculate the finite temperature thermal expan- 
sion and heat capacity, a quasi harmonic approximation 
(QHA) was employed [H [ll|. Using the DFT resuhs 
for the phonon dispersion relation and minimizing the 
QHA free energy with respect to the lattice parameter 
{^\t = 0), they found the linear thermal expansion of 
the lattice parameter -a- for graphite, GE and diamond. 
The linear thermal expansion for the lattice parameter 
of GE, 7a = where a — \/3ao, was found to be 

negative and about 7a — —0.35 x 10^^ K^^ at room 
temperature. QHA+GGA includes only weakly anhar- 
monicity but it includes quantum effects and zero point 
energy. The strong anharmonic coupling of the bending 
and stretching modes in GE, which are essential for the 
GE stability, are accurately described by classical em- 
pirical potentials at non-zero temperatures [ll]. In our 
study the negative thermal contraction is reported for the 
length -L- of GA and GE samples. Since ripples appear 
in GE and GA membranes, we expect a contraction of 
the length of the systems. 7l is not directly related to 
7a. In fact, dL is the variation of the measured length 
along x-direction while da is the variation of the lattice 
parameter not necessarily in the x — y plane. Therefore 
our results for the contraction of the GE's length are not 
directly related to those of QHA-I-GGA for contraction 
of the lattice parameter. 

In a recent experiment (l7| the morphology of differ- 
ent graphene membranes (micron size membranes) sus- 
pended across trenches on Si/Si02 substrates were in- 
vestigated using AFM and STM. TECs were found to 
be about ^exp ~ —0.7 x 10~^ at room tempera- 
ture. With increasing temperature the measured TECs 
approaches zero nonlinearly. Here, we find a similar 
contraction in the longitudinal size (e.g. 7l ~ —9.0 x 
IQ-^R-^ at T=300K for GE) which approaches zero 
with increasing temperature but which is an order of 
magnitude larger. To compute the thermal contraction 
we calculated the longitudinal size (not the arc length as 
done in the experiment flT'l). There are important dif- 
ferences between our approach and those from the above 
mentioned experiment. First, we use the initial length (as 
a reference for the length) of our samples as the length of 
the sheet at zero temperature. Second, the sample size is 
different. Third, there is no substrate in our study which 
could reduce the contraction. Our suspended sample is 
free to drift in the x—y plane while in the experiment, the 
substrate prevents free drift in the x — y plane in addition 
to movement along the z-direction which causes a lower 
contraction. Notice that < lexp < la- We found that 
in thermal equilibrium and under supported boundary 
conditions, the equilibrium size of GE is longer than GA. 
Therefore, the GA surface is much more corrugated than 
GE. 

The larger contraction of GA is due to the larger am- 



plitude of the ripples as compared to GE. The reason 
is that hydrogen atoms which are below and above the 
sheet attract or repel the carbon atoms and push them 
in different directions, (i.e. random thermal fluctuations 
of the hydrogens at finite temperature). Therefore, we 
expect an increasing randomness of the GA sheet and 
expect the formation of different patterns of ripples in 
a GA sheet and a larger corrugation. The thickness 
of GA is larger than GE, hence one expects for such a 
thicker material, ordinary positive thermal expansion (or 
at least smaller thermal contraction with respect to GE) . 
However, we found negative thermal contraction for GA 
which is larger than the one for GE. It is interesting to 
note that above 1500 K we observe a buckling of GA. Fig- 
ure |3] shows a snap shot of a buckled GA sheet at 2000 K. 
A similar buckling was observed in experiments on sus- 
pended GE 1171 and in a circular GE sheet subjected to 
radial strain [l8|. Therefore, increasing temperature in- 
duces larger axial strains in GA. 



B. Roughness 

Figure [2Kb) shows the average square root out of plane 
deviation of the carbon atoms (i.e. static roughness, i.e. 
w = \J {K^) — {h)"^) versus time for three typical tem- 
peratures. For long times w fluctuates around 0.2 A for 
GE and 1.3 A for GA which is 6.5 times larger. Fig- 
ure [JKc) shows the variation of (/i^) (where average is 
over lateral size of the system) and (/i^) (where aver- 
age is over longitudinal size of the system) for GE and 
GA at T=700K. Because of the specific boundary con- 
dition the ripples appear in the cc-direction, while in the 
y-direction the fluctuations are much smaller and for GE 
they are almost zero. This leads to a larger randomness 
in GA which relates to a less stiff material as compared 
to GE Q. This is illustrated in Fig. 2] which shows the 
corresponding contour plot for GA and GE at T=300K 
and 700 K. Notice the difference in scale for the ampli- 
tude. For example for GA at 700 K the amplitude of the 
ripples are in the [-4,4] A range and [-2,2] A for GE. In 
the lateral boundaries we see larger amplitude variations. 
The free boundaries in GE exhibit larger free vibrations 
than GA. This can also be inferred from the top panel 
of Fig. [2Kc). Notice that the typical wave lengths (A) for 
the ripples in GA are very different as compared to the 
one of GE. Figure [5] shows the amplitude of the Fourier 
transform oi hx, {\hq\'^) (with an ensemble average taken 
over 60 samples) , for GA and GE as function of the wave 
length along the a;-direction (taken in the middle of the 
system around y—0), Qx — ^ for two different temper- 
atures. Independent of temperature, the ripples in GA 
span a wider range of wave lengths, for example they span 
length scales in the (30-125) A and (25-140) A ranges, for 
T = 300K and T = 700 K, respectively. This compares 
with ripples in GE having typical length scale around 
60 A for T = 700K while for T = 300K the membrane is 
almost flat. 
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FIG. 2: (Color online) (a) Temperature dependence of the thermal contraction coefficient of a square sheet of GA and GE 
subjected to longitudinal supported boundary and lateral free boundary conditions. Solid curves are fits according to a + P/T. 
(b) Roughness of GA and GE as a function of time, (c) The variation of (h^) (averaged over lateral size) and {hy) (averaged 
over longitudinal size) at 700 K. 
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FIG. 3: (Color online) Buckling of graphane at 2000 K. 
Dashed curve shows the convex shape of graphane. 
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FIG. 4: (Color online) Contour plot of the z-position of the 
carbon atoms of GA and GE at 300 K and 700 K correspond- 
ing to the situation of Fig. EJc). The lateral edges are free 
while the longitudinal edges are not allowed to move in the 
z-direction. 



FIG. 5: (Color online) The amplitude of the Fourier trans- 
form of hx of GA and GE atoms along the a;-direction (av- 
eraged over 60 samples) for two temperatures T = 300 K and 
T = 700K (the inset). 



Next we estimate the roughness exponent x of GA 
(where only the carbon atoms are considered) which is 
obtained from the second order structure function, i.e., 
S{d) = \ {h{x+6)-h{x))\'^ which scales as [H where S 
counts the atoms along the a;-direction, e.g. (5 = 4 refers 
to the fourth atom in the armchair direction {S is typi- 
cally 2 A). Figure [S^a) shows the variation of S{S) versus 
S for GA and GE for two temperatures. Notice that S{d) 
for GA is almost an order of magnitude larger than for 
GE indicating the larger corrugation. S{S) increases with 
S up to some critical 5- values, Sc, which is related to the 
above mentioned characteristic length scale of the rip- 
ples. For 6 < 5c the height of the atoms are correlated 
while they become uncorrelated for S > Sc- The slope 
of S{S) gives X whose temperature dependence is shown 
in Fig. E^b). Below room temperature GE has a larger 
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FIG. 6: (Color online) (a) Variation of the structure function 
versus S (log-scale) for GE and GA at 300 K and 700 K. Here 
S counts the atoms along armchair direction and is about 
2 A. Dashed lines are fits to 5^. (b) Roughness exponent as 
a function of temperature for GA and GE. 
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FIG. 7: (Color online) (a) Variation of the potential energy 
(Ep) per atom versus temperature for both C atoms in GA 
and GE and H atoms in GA. The right scale is for H atoms 
and the left scale is for C atoms, (b) Variation of total energy 
(Et = Ep + Ek) per atom versus temperature for both GA 
and GE. The right scale is for GA. The error bars are less 
than 0.005 eV/atom. 



X (~1.2) which implies that in this temperature range 
the GE membrane has a smoother surface as compared 
to GA. However at high temperatures both systems ap- 
proach the situation with random height fluctuations, i.e. 
X ^ 1. Indeed, one expects that the presence of the sp^ 
bonds in GA (making GA effectively a much thicker ma- 
terial than GE) decreases the roughness with respect to 
GE while we found it increases. 



C. Heat capacity 

Figure [TJa) shows the variation of the potential energy 
per atom {Ep) for C atoms in GE, C atoms in GA and 
H atoms in GA versus temperature (all data points were 
obtained by averaging over an ensemble of 150 samples). 
The potential energy is a measure of the binding energy 
(or bond energy), i.e. BE — —Ep. The binding en- 
ergy of the C atoms in GA is larger than the one in GE. 



The right scale in Fig. El^a) is for the potential energy 
of H atoms in GA. The variation of the potential energy 
of H atoms with temperature is smaller than the one 
for C atoms in GA. Note that an accurate study of the 
GA binding energies at low temperature needs ab-initio 
molecular dynamic simulations where quantum fluctua- 
tions are included. At low temperature, where the dif- 
ference between energy of two quantum states of the sys- 
tem is larger than the thermal energy, many body effects 
and the chemical covalent bond energies become sensi- 
tive to the quantum states. Here we found that the sum 
of the binding energy of C-H bonds at 50 K is 5.012 eV 
per atom which is comparable to the recently reported 
value at T=OK, i.e. 5.19 eV per atom using DFT cal- 
culations [1^. The extrapolation of our data to r=0K 
gives 5.02 eV per atom which indicates that quantum cor- 
rections are of order 0.17eV per atom or about 3.3%. 
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Extrapolation in Fig. ^a), for the binding energy of C 
atoms in GE gives 7.33eV/atoni (for all reported energies 
the error bars are less than 0.005eV/atom). 

Figure ^h) shows the total energy per atom (the 
sum of potential and kinetic energy of C and H atoms 
Et — Ep + Ek) versus temperature (potential energies 
are also shown in Fig. [7Kb) for comparative purposes. 
The right scale in Fig.lTJb) is for the total energy of GA. 
At a temperature around 1500 K we observe evaporation 
of H atoms at the free edges. The total energy varies 
linearly with temperature and gives the corresponding 
lattice contribution to the molar heat capacity at con- 
stant volume (the average size of the system after re- 
laxation is taken constant) Cv — ^^§f- which for GE 
is 24.98±0.14J/molK and for GA is 29.32±0.23 J/molK. 
This is comparable to the proposed classical heat capac- 
ity at constant volume i.e. Cy = 35R ~ 24.94 J/molK, 
i.e. the Dulong-Petit limit where 5R is the universal gas 
constant. The larger heat capacity for GA is due to the 
extra storage of vibrational energy in the C-H bonds. Our 
result for GE is in agreement with those obtained from 
Monte Carlo simulations of Ref. _21!] but with smaller er- 
ror bars. The reported heat capacity for GE in Ref. [2lj 
is Cv ^ 25±1.0 J/molK for r=1000K. The used size for 
the GE samples in Ref. [2l| are 13% smaller than those 



in the present study. 



IV. CONCLUSION 

In this paper we studied the thermal properties of sus- 
pended graphane and compared the obtained results to 
the one found for graphene. We found that the rough- 
ness, amplitude and wave lengths of the ripples are very 
different. The thermal contraction effect for graphane 
is larger than for graphene. Above 1500 K we found 
that graphane is buckled and starts to loose H-atoms 
at the edges of the membrane. Roughness in graphane is 
larger than in graphene and the roughness exponent in 
graphene decreases versus temperature (from 1.2 to 1.0) 
while for graphane it stays around 1.0 implying random 
uncorrelated roughness. Fourier analysis of the height of 
the C atoms showed that the ripples in graphane exhibit 
a larger range of length scales as compared to the one 
for graphene. Heat capacity of graphane is found to be 
14.8% larger than the one for graphene. 
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